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ABSTRACT 

We present results of a survey for giant Lya nebulae (LABs) at z = 3 with 
Subaru/Suprime-Cam. We obtained Lya imaging at z = 3.09±0.03 around the SSA22 
protocluster and in several blank fields. The total survey area is 2.1 square degrees, 
corresponding to a comoving volume of 1.6 x 10 6 Mpc 3 . Using a uniform detection 
threshold of 1.4 x 10~ 18 erg s _1 cm~ 2 arcsec~ 2 for the Lya images, we construct a 
sample of 14 LAB candidates with major-axis diameters larger than 100 kpc, including 
five previously known blobs and two known quasars. This survey triples the number 
of known LABs over 100 kpc. The giant LAB sample shows a possible "morphology- 
density relation" : filamentary LABs reside in average density environments as derived 
from compact Lya emitters, while circular LABs reside in both average density and 
overdense environments. Although it is hard to examine the formation mechanisms of 
LABs only from the Lya morphologies, more filamentary LABs may relate to cold gas 
accretion from the surrounding inter-galactic medium (IGM) and more circular LABs 
may relate to large-scale gas outflows, which are driven by intense starbursts and/or 
by AGN activities. Our survey highlights the potential usefulness of giant LABs to 
investigate the interactions between galaxies and the surrounding IGM from the field 
to overdense environments at high-redshift. 

Key words: galaxies: formation - galaxies: evolution - cosmology: observations - 
cosmology: early universe 



1 INTRODUCTION 

Lya blobs (LABs) are spatially extended Lya 

nebulae seen in the high - redsh if t Universe (e.g. , 

19961: iKeel et al ' " 



Francis et al 



1999; 

■rr 



Steidcl et al 



2000; 
2005; 



Palunas et all 12004: Matsuda et al. 2004; Dev et al 
Nilsson et al.ll2006l:lGreve et al.ll2007l:ISmith fc Jarvisll2007l : 
Prescott. Dev. fc Jannuzil 120091 : iVang et all 120091 ). LABs 



are thought to relate to the formatio n of massive galaxies 
l|Dev et al. I l2005l : iMatsuda et all 120061 ') and to be indicative 
of strong interactions between the inter-galactic medium 



* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan, 
f E-mail: yuichi.matsuda@durham.ac.uk 
| Carnegie fellow. 



(IGM) and gala xies with intense star-f ormation activities 
and /or AGNs jFurlanetto et alj l2005l b To explain the 
formation mechanisms of LABs, at least three possible 
ideas have been proposed: cold gas accretion, galactic 



dHaiman, Spaans, & C 


uataert 


l200d; Chapman et al. 


200 lib 



Taniguchi fc Shioval 
In spite of extensive ob- 
servational and theoretical efforts in the decade after 
the first discovery of LABs, th e formation mechanisms 
of LABs ar e still controversial |Mori fc Umemural 120061: 
i Geach et alj 120091: iDiikstra fc Loebl 120091: iGoerdt et al 



bold : iFaucher-Giguere et all l20ld : Isii imizu fc Umemural 

bold) . 

Among the LABs, special attention have been given to 
the largest examples with the spatial extents of ~ 100 — 
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Table 1. Summary of narrow-band observations 



Field 


RA 


Dec 


Date 


Exposure 


Area 


FWHM 




Depth 




(J2000) 


(J2000) 


(mm/yyyy) 


(hours) 


(arcmin 2 ) 


(arcsec) 


(cgs) a 


(ABmag) 6 


SXDS-C 


02:18:00.0 


-05:00:00 


08, 09, 10/2005 


5.2 


682 


1.0 


0.81 


26.3 


SXDS-N 


02:18:00.0 


-04:35:00 


10/2005 


4.8 


740 


1.0 


0.94 


26.2 


SXDS-S 


02:18:00.0 


-05:25:00 


08, 10/2005 


4.8 


737 


1.0 


0.82 


26.3 


GOODS-N 


12:37:23.6 


+62:11:31 


04/2005 


10.0 


869 


1.1 


0.69 


26.6 


SDF 


13:24:39.0 


+27:29:26 


04/2004, 04/2005 


7.2 


805 


1.0 


0.67 


26.5 


SSA22-Sbl 


22:17:34.0 


+00:17:01 


09/2002 


7.2 


633 


1.0 


0.92 


26.3 


SSA22-Sb2 


22:16:36.7 


+00:36:52 


08/2004 


5.5 


487 


1.0 


0.96 


26.3 


SSA22-Sb3 


22:18:36.3 


+00:36:52 


08, 09/2005 


5.5 


537 


1.0 


0.89 


26.3 


SSA22-Sb4 


22:19:40.0 


+00:17:00 


08, 09, 10/2005 


5.5 


529 


1.1 


1.15 


25.9 


SSA22-Sb5 


22:15:28.0 


+00:17:00 


09/2005 


5.5 


565 


1.0 


1.06 


26.1 


SSA22-Sb6 


22:14:30.7 


+00:33:52 


10/2005 


5.5 


572 


1.0 


0.92 


26.3 


SSA22-Sb7 


22:17:42.7 


+00:56:52 


09, 10/2005 


5.5 


480 


1.0 


1.02 


26.2 



"The 1-<t surface brightness limit (10 18 erg s 1 cm 2 arcsec 2 ). 

6 The 5-<r limiting magnitude calculated with 2 arcsec diameter aperture photometry. 



200 kpc (hereafter giant LABs) because of their spectac- 
ular morphologies and possible association with protoclus- 
ters dSteidel et al.ll200d: IPalunas et al.ll2004l : IPrescott et~al] 
feOOSUMatsuda et al.ll2009h . At present, there are only six 
known giant LABs over 100 kpc, and they have been se- 
lecte d by using different quality data sets and different meth- 
ods llFrancis et al.lll996l: ISteidel et al.ll2000l : iDev et al.ll2005l : 
iGreve et al.ll2007l : llViatsuda et al.ll2009r ). It is therefore diffi- 
cult to examine their statistical properties. In order to con- 
struct a statistically reliable sample of giant LABs and so 
test their possible association with overdense environments, 
we undertook a deep, wide-field Lya imaging at z = 3.1. 

In this letter, we use AB magnitudes and adopt cosmo- 
logical parameters, JIm = 0.3, = 0.7 and Ho = 70 km 
s _1 Mpc -1 . In this cosmology, the Universe at z = 3.1 is 2.0 
Gyr old and 1.0 arcsec corresponds to a physical length of 
7.6 kpc at z = 3.1. 



2 OBSERVATIONS AND DATA REDUCTION 

We briefly describe our observations and data reduction, 
although the details will be reported in a separate paper 
(Yamada et al. in prep). The summary of the observa- 
tions and data is listed in Table 1. The imaging observa- 
tions were carried out betwee n September 2002 and Octo- 
ber 2005 using Suprime- Cam dMiyazaki et all 120021 ') on the 
8.2-m Subaru Telescope l|lve et al.ll2004l ). Suprime-Cam has 
a pixel scale of 0.2 arcsec and a field of view of 34 x 27 
arcmin 2 . We obtained narrow-band (NB497) images for 
12 pointings: Great Observatories Origins Deep Survey- 
North (GOODS-N), Subaru Deep Field (SDF), three fields 
in Subaru-XMM Deep Survey (SXDS-C, N, and S) and seven 
fields around SSA22a (SSA22-SM-7). The SSA22-SM field 
was the first field of our survey and centred at SSA22a, 
which cont ains the prot o cluste r region at z = 3.09 dis- 
covered by ISteidel et all (|2000T l. Initial results of the ob- 
servations in the SSA22-Sbl have been alre a dy publishe d 
|Havashino et all 120041 : iMatsuda et all |2004 120051 . 120061 ). 
The NB497 filter has a central wavelength of 4977 A and 
FWHM of 77 A which corres ponds to the redshift range 
for Lya at z = 3.062 - 3.126 (jHavashino et al.1 120041 '). The 
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Figure 1. Colour magnitude plot of BV - NB497 vs. NB C 
for NB C detected sources (black dots). The solid line represents 
colour criterion (BV — NB497 = 0.7) used for narrow-band excess 
objects (larger dots). The blue squares and red circles indicate 
giant (major-axis diameters^ 100 kpc) LAB candidates without 
QSO and with QSO, respectively. All magnitudes and colours are 
measured with isophotal apertures. 



width of the redshift slice is 59 comoving Mpc. For the 
SSA22 fields, we obtained broad-band (B and V) images 
in our observing runs. For GOODS-N, we used archival raw 
B and V-band images l|Capak et alj|2004l ). For the SDF and 
SXDS fields, we used public, reduced B and I /-band images 
ijKashikawa et al.ll2004l : iFurusawa et al.ll2008h . 

We reduced the ra w data with sdfred (jYaei et al.1 
120021 ; lOuchi et all 120031 ) and iraf. We calibrated the as- 
trometry of the i mages using the 2 MASS All-Sky Catalog 
of Point Sources ijCutri et al.ll2003r ). For photometric cal- 
ibration, we used the photometric and spectrophotomet- 
ric standard stars, SA113, SA115, FEIGE34, Hz44, P1 77D 
GD248, SA95-42, L DS749B, BD+332642, and G24-9 ijOkd 
ll99d;lLandoitlll992l '). We corr ected the magnitudes using the 
Galac tic extinction map of ISchlegel. Finkbeiner. fc Davis! 
(1998). We aligned the combined images and smooth with 
Gaussian kernels to match their seeing to a FWHM of 1.0 
or 1.1 arcsec depending on the original seeing. We made 
BV images [BV = (2B + V^)/3] for the continuum at the 
same effective wavelength as NB497 and made NB C (con- 
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Table 2. Properties of the 14 giant LAB candidates 



ID 


RA 


Dec 


a a 


Area 


^Lyo 


F b 


<Slae 


^■spec 


Note 




(J2000) 


(J2000) 


(kpc) 


(arcsec 2 ) 


(10 43 erg s" 1 ) 










SSA99 Shi T,AR1 


99-1 7-9^ 


-1-00-1 9-^7 7 


175 


181 + 14 


8.1 + 0.6 


56 


2.7 


3.099 c 


Si / 1 TYl ^ / £2 1 1 V* rvi TY1 3 
OLil.il j 3 LI U111II1 


SSA22-Sb6-LAB1 


22:13:48.30 


+00:31:32.8 


166 


116 + 9 


5.8 + 0.4 


0.69 


0.6 


3.094 d 




SSA22-Sbl-LAB2 


22:17:38.99 


+00:13:27.8 


157 


137 + 8 


6.8 + 0.3 


0.59 


3.7 


3.091 c 


X-ray fe /8^m l 


SSA22-Sb5-LAB1 


22:15:33.56 


+00:25:16.9 


147 


59 + 7 


3.8 + 0.4 


0.80 


-0.5 






SSA22-Sb3-LAB1 


22:17:59.45 


+00:30:55.7 


126 


102 + 8 


20.4 + 0.3 


0.52 


1.2 


3.099 e 


QSOVRadio' 


GOODS-N-LAB1 


12:35:57.54 


+62:10:24.9 


124 


47 + 7 


5.4 + 0.5 


0.77 


0.9 


3.075^ 


QSO^'/X-ray" 1 


SSA22-Sb2-LAB1 


22:16:58.37 


+00:34:32.0 


121 


60 ± 15 


2.0 + 0.6 


0.70 


1.2 






SSA22-Sb2-LAB2 


22:16:56.40 


+00:27:53.3 


115 


48 ± 11 


1.4 + 0.2 


0.73 


-0.1 






SSA22-Sbl-LAB5 


22:17:11.66 


+00:16:44.4 


110 


43 ± 11 


1.3 + 0.3 


0.74 


1.0 




8/im 1 / submm" 


SSA22-Sb5-LAB2 


22:15:30.27 


+00:27:43.6 


107 


53 + 7 


2.1 + 0.3 


0.66 


-0.1 






SSA22-Sb6-LAB4 


22:14:09.58 


+00:40:54.6 


107 


32 + 4 


2.0 + 0.2 


0.79 


-0.1 


3.116 d 




SSA22-Sbl-LAB3 


22:17:59.14 


+00:15:28.7 


103 


75 + 9 


5.2 + 0.2 


0.48 


1.7 


3.096 9 


X-ray ! 


SXDS-N-LAB1 


02:18:21.31 


-04:42:33.1 


101 


68 + 5 


3.3 + 0.2 


0.51 


-0.4 






SSA22-Sbl-LAB16 


22:17:29.01 


+00:07:50.2 


101 


28 + 8 


0.8 + 0.2 


0.80 


-0.2 


3.104'' 


X-ray ! /8^tm l /submm n 



a Maj or-axis diam e ter, b Fi lamentarity (F = fo r a circle, F = 1 for a fi lament, see tex t for m ore detaill. ISteidel et al 1 ll2003h. this 
work ' IShen et all J2007I1 . ^Bareer et ail ^0021), ^Matsuda et al.l ||2005|1, 'jMatsuda et alj I2006T), iWebb et alj (1200911 . -^Chapman et all 



l|200ll l. * lBasu-Zvch fc Scharj l|2004h . ICondon et al.l l|l998l l. " lAlexander et alj J2003h . iGeach et alj J2005I '). jGeach et alj J200a ) 






Figure 2. Pseudo-colour images (B for blue, NB497 for green, V for red) of the 14 giant LABs. The size of the images is 40 X 40 arcsec 2 
(~ 300 X 300 kpc 2 ). The yellow contours indicate isophotal apertures with a threshold of 1.4 X 10 -18 erg s — 1 cm -2 arcsec -2 . The white 
horizontal bar in the lower right image represents the angular scale of 100 kpc (physical scale) at z = 3.1. 



tiriuum subtracted NB497) images for emission-line images. 
The total survey area after masking low S/N regions and 
bright stars is 2.12 square degrees and the survey volume is 
1.6 x 10 6 com oving Mpc 3 . Th i s is 12 times larger than the 



survey area of iMatsuda et al l d2004l ) and 100 times larger 
than that of ISteidel et all (|200ol )~The 1-a surface bright- 



ness limits of the NB C images are 0.7 — 1.2 x 10 



" 1S erg s -1 



3 RESULTS 

Object detection and photometry are performed using 
the double image mode of SExtractor version 2.5.0 
( Bcrt in fc Arnoutsl [l996). For source detection, we use 
smoothed NB C images with Gaussian kernels to match their 
seeing to a FWHM of 1.4 arcsec in order to slightly increase 
the sensitivities for diffuse extended sources and to make 
all the images the same seeing size. We use the same detec- 
tion threshold (DETECT-THRESH) of 1.4 x 10 -18 erg s -1 
cm -2 arcsec -2 (or 28.5 mag arcsec -2 ) for all the 12 fields. 



The magnitudes and colours are measured with isophotal 
apertures defined in the NB C images. 

In Fig. 1, we plot the BV - NB4Q7 colours and NB C 
magnitudes of the A^B c -detected sources. The solid line rep- 
resents the colour criterion used for narrow-band excess ob- 
jects, BV — NB4Q7 — 0.7, which corresponds to an observed 
equivalent width of EW h s = 80 A. From these narrow-band 
excess objects, we make a diameter-limited catalog of 14 
LABs down to a major-axis diameter of the isophotal aper- 
ture of a ^ 13 arcsec (or ^ 100 kpc at z = 3.1). 

For the LAB selection, we use the major-axis diam- 
eters rather than isophotal area, in order to cover LABs 
with asymmetric structures. For example, cold stream mod- 
els pre dicted that LABs have asymmetric, long and thin fila- 
ments (|Goerdt et al.ll20ld ; lFaucher-Giguere et alj|20ich . An 
alternative quantity for the LAB selection may be Lya lu- 
minosity. However, the Lya luminosity could be dominated 
by a bright central core, such as starbursts in the central 
galaxy and AGN. 

Six out of the 14 LABs have been spectroscopically con- 
firmed by previous surveys l|Steidel et al.|[200ol : iBarger et al.l 
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Figure 3. Sky distribution of the 14 giant LABs and smoothed dens ity maps of ~ 2000 compact LAEs at z ~ 3.09. I n the left panel 
(a), th e small black box indicates SSA22a field bv lSteidel et al] d200Cl . S00) and the dashed box indicates SSA22-Sbl by Matsuda et all 
|2004 M04). The thick bars show the angular scale of 20 comoving Mpc at z = 3.1. The blue squares and red circles indicate the giant 
LABs without QSO and with QSO, respectively. The contours represent LAE overdensity, <5lae = (n — n)/n =0, 1, 2, 3, 4, 5, and 6. 



2002; lMatsuda et al. 2005, 120061 ; IShen et alj |2007). For two 
new LABs, we carried out spectroscopic follow-up observa- 
tions with Magellan/IMACS in June 2009, and identified 
Lya emission from these LABs (see Ono et al. in prep for 
more details). 

The properties of the 14 giant LABs are listed in Ta- 
ble 2. We rename the LABs in the new sample since ini- 
tial surveys. SSA22- Sbl-LAB16 was named as LAB 18 in 
iMatsuda. et~ai1 fa004) . Note that SSA22-SM-LAB1 discov- 
ered by ISteidel et al.l |2000r i is still the largest one in this 
sample. SSA22-Sb3- LAB1, and GOODS-N-LAB1 are as - 
sociated with QSOs jBarger et alj|2002l ; IShen et all 12007ft . 
We apply the "QSO" labels only to optically bright known 
QSOs but some of the other LABs have potential signs 
of o bscured AGN. All the 5 LABs from t he initial sur- 
veys (|Steidel et al.ll200d ; lMatsuda et al.ll2004T ) were detected 
at X-ray, 8/im, and/or submm wavelength follow-up ob- 
servations (IChapman et al. 200ll; Basu-Zvch fc Scharj|2004 
iGeach et al.ll2005l . 120091 ; IWebb et al.ll2009l V Thus, although 
there are only two LABs are apparently associated with 
QSOs, the AGN fraction of the 14 LABs may be higher 
when accounting for more obscured AGN. 

The thumbnail images (40 x 40 arcsec 2 ) of the 14 giant 
LABs are displayed in Fig. 2. The LABs show a wide variety 
of Lya morphologies. While some LABs appear to have cir- 
cular shapes, some have filamentary (or elongated) shapes. 
We quantify the Lya morphology by defining "filamentar- 
ity", 



F 



((isophotal area)/(7r x (a/2) )) 



where a is the major-axis diameter. For example, a circle 
has F = and an extremely thin filament has F = 1. The 
nlamentarities of the LABs range F ~ 0.4 — 0.8. We esti- 
mate the uncertainties of the Lya properties by putting the 
thumbnail NB C images of each LAB at 100 random positions 
on the original NB C images and measuring the deviations. 

Fig. 3 shows the spatial distributions of the 14 giant 
LABs and smoothed density maps of ~ 2000 Lya emitters 




Figure 4. Filamentarity of the 14 giant LABs as a function of 
the overdensity of LAEs. The blue squares and red circles indicate 
giant LABs without QSO and with QSO, respectively. The error 
bars show 1-cr uncertainties. The filamentarity of the LABs shows 
a weak anti-correlation with the overdensity of LAEs. 



(LAEs) selected by Yamada et al. (in prep) with the same 
data. Since some bright Lya knots in the giant LABs are 
also detected as single or multiple LAEs, we exclude such 
LAEs from the LAE sample in this analysis. We make the 
density maps by smoothing the LAE spatial distributions 
with a Gaussian kernel of a = 1.2 arcmin, or FWHM = 5.3 
comoving Mpc. The smoothing kernel size is chosen to match 
the median distance between the nearest neighbours in the 
LAE samples in the blank fields (SXDS, GOODS-N, and 
SDF). The contours represent LAE overdensities, 5lae = 
(n — n)/n =0, 1, 2, 3, 4, 5, and 6, where n is the mean LAE 
number density in the blank fields. 

In Fig. 4, we plot the nlamentarities of the 14 LABs as a 
function of the LAE overdensities. The filamentarity shows 
a weak anti-correlation with the LAE overdensity. While 
more filamentary LABs reside in average density environ- 
ments, more circular LABs reside in both average density 
and overdense environments. The Spearman's rank corre- 
lation coefficient is r s — —0.56. We can rule out random 
distributions with 96% confidence. We have confirmed that 
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the results do not change significantly if we use isophotal 
area for the LAB selection, suggesting that the correlation 
is not due to the selection method. 



4 DISCUSSION AND CONCLUSIONS 

Based on deep, wide-field Lya imaging, we construct a sam- 
ple of 14 giant LAB candidates at z — 3.1 from a volume of 
1.6 x 10 6 comoving Mpc 3 . This is the largest sample of giant 
LABs and triples the number of known LABs over 100 kpc. 
Our giant LAB sample shows a wide variety of Lya mor- 
phologies and resides not only in overdense environments, as 
derived from LAEs, but also in lower dense environments. 
We find a possible hint for "morphology-density" relation of 
the LABs: the Lya filamentarity seems to differ as a function 
of the local density environments. 

How can we interpret this possible morphology-density 
relation of the LABs? The Lya morphology may relate to 
the formation mechanisms of LABs. According to recent nu- 
merical simulations, more filamentary LABs may be good 
candidates for cold gas accretion from the surrou nding IGM 
(|Goerdt et al.ll20ld ; [Faucher-Giguere et aUl201Cft . Although 
direct evidence for such gas inflows is not foun d around star- 
forming galaxies at z ~ 2 (|Steidel et al.ll20lbT ). recent stud- 
ies of the metallicity of star-forming galaxies from low- to 
high-redshifts indicate that gas inflows m ay still be dom- 
inant in the field environment at z > 3 (tMannucci et al.l 
2010). More circular LABs may relate to large-scale gas out- 
flows, whi ch are driven by intense starbursts and/or AGN 
activities (M on fc Umemurall2006h . At high-redshift, star- 
formation and AGN activities in overdense environments are 
known to be se veral times higher than those in the field envi- 
ronments fe.g.. lSmail et aLll2003T ). Future spectroscopic and 
multi-wavelength follow-up observations would enable us to 
investigate the gas dynamics and the variations of the star- 
formation and AGN activities in giant LABs as a function 
of the environments and to test the interpretations. 
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